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PRELIMINARY EXPERIMENTS ON FORCES AND MOMENTS OF AN 
OSCILLATING WING AT HIGH-SUBSONIC SPEEDS 
By S. A. Clevenson and E. Widmayer, Jr. 

SUMMARY 


Experimental results of lifts and moments about the quarter chord 
of an oscillating wing at high— subsonic Mach numbers are presented. A 
comparison of the experimental magnitude of the lift vector with the 
theory as given by Dietze showed good agreement. Comparisons of the 
moments and the quadrature component of lift with theory indicated that 
some refinements in the testing technique are necessary for the experi- 
mental determination of these quantities in the transonic range. 


INTRODUCTION 


An extensive amount of work has been done on the theory of air 
forces on a harmonically vibrating wing in a two-dimensional compres- 
sible medium by Possio, Frazer, Dietze (references 1 to 3), and others. 
Flutter of wings in a compressible medium has also been treated by 
Garrick (reference L), Frazer and Skan (reference 5)> and others in 
which the aforementioned theoretical air forces are used. Numerous 
wind— - tunnel tests have also been made on the flutter of cantilever wings 
in the subsonic— speed range (reference 6 ). These tests 3 hov reasonable 
agreement with the indicated theoretical trends. Some work on the 
experimental determination of the oscillating air forces at high- 
subsonic speeds has been done by Bratt and Chinneck (reference 7 ), and 
the oscillating moments on a small airfoil oscillating about the mid- 
chord position were determined at high— subsonic speeds for a range of 
reduced frequencies of 0.0025 to 0.015. The aerodynamic coefficients 
have been calculated by Dietze for a reduced frequency range to = 0.6 

and for a Mach number range to M = 0 . 7 . . 
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The present paper deals with some experiments at high— subsonic 
speeds in which a two-dimensional wing was oscillated about the quarter- 
chord position in a reduced frequency range from 0.15 to 0.35. Since 
extension of the theory into the transonic range is extremely difficult 
and open to considerable doubt because of the occurrence of shocks and 
mixed flow on the wing, it appears desirable to perform an experiment 
for determining the oscillating air forces and moments in this region. 
Although this paper presents few data in the transonic region, it does 
present an experimental method for obtaining these forces and moments. 

A comparison is made of the experimental forces and moments with the 
theory of Dietze. Although comparisons would not be made for all 
components, it is believed that the presented method can be refined to 
give sufficient accuracy for the determination of the components of lift 
and moment of an oscillating airfoil in the transonic range. 


SYMBOLS 


b 

l 

P 

V 

Z 

L 

Q 

M 

a 

% 


semichord 

representative length (2b) 

mass density. of gas 

velocity of gas 

unit span of wing 

oscillating lift (L Q sin(o>t + <P]_)) 

oscillating moment about quarter chord (Q 0 sin(o>t f cpg)) 
Mach number 

amplitude of pitching oscillation (<Xq sin cut) 
natural circular bending frequency 


CD 

da 


natural circular torsion frequency 

circular frequency of oscillation 
reduced frequency (bcu/v) 

slope of lift curve 
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<Pl 

$2 

K* 


SD 


K" 


SD 


K' 


DD 


phase angle "between torsional displacement and lift 
phase angle "between torsional displacement and moment 
lift coefficient in phase with torsional displacement 

r 

L 0 cos Cf>]_ CD^ 
v rtpV 2 hZa 0 2 

lift coefficient in quadrature to torsional displacement 
f L 0 sin cpj N 

rr p V 2 "bZa 0 

moment coefficient about quarter chord in phase with torsional 


[ Q n COS (Do o pi 

displacement | — — + £ orT { 

itpV 2 b 2 Za ° 


K" 


DD 


moment coefficient about quarter chord in quadrature with 

t \ 

Q 0 sin qp g 


torsional displacement 


k *pV 2 b 2 Za 


The subscript o represents the vector magnitude. 


APPARATUS 

Model 


The wing has a span of 2k inches and a chord of 8 inches, with an 
NACA 65 —OIO airfoil section. It was of fabricated construction, having 
a steel spar with duralumin skin and was mass balanced about the 
quarter chord to eliminate any mass— inertia force. Thi 3 airfoil was 
supported on bearings at each end and had only freedom to pitch about 
the quarter-chord point. The wing oscillated about a mean angle of 
attack of 0° and had a maximum amplitude of ±2.30°. Two variable- 
frequency, electromagnetic shakers phased to produce pure moment and 
a torsion rod were attached near one end of the wing (fig. 1). The 
other end of the torsion rod was fastened to the tunnel wall. Various 
resonant frequencies could be obtained by the use of torsion rods 
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having different torsional stiffness. .The frequencies of the wing— rod 
combinations used are as follows: 


Wing— rod combination 

% 

( radians / sec ) 

01- 

a 

(radians/sec) 

1 

528 

194 

2 (Data not reported) 

546 

275 

3 (Data not reported) 

572 

361 


Tunnel 


The test section' of the Langley 4|— foot flutter research tunnel 


was modified by means of auxiliary walls so that the test section was 


24 inches wide. 


A 


_1_ 

16 


—inch gap existed between the model and the tunnel 


walls. All tests reported were made in Freon— 12 in which the speed of 
sound was about 520 feet per second. The use of this test medium 
permits the attainment of more than twice the reduced frequency as 
obtained in air for a given Mach number and torsional frequency. 

Reynolds numbers for these tests were from 1 x 10^ to 5.5 x 10^. 


Instrumentation 

The instrumentation may be divided into two groups: the first’ 

indicating the motion of the wing in the tunnel and its reaction to the 
flowing gas and the second indicating the supplied power necessary to 
maintain the oscillation. A 14— channel recording oscillograph was used 
to record and correlate data from both groups with respect to time. 

The -time history of the pitching motion of the wing was obtained 
from the signal of an electric strain gage which measured the torsional 
stress in the torque rod. A light beam, reflected from a mirror 
fastened to the end of the wing and focused on a translucent scale, also 
gave an indication of the torsional amplitude of oscillation (fig. 1). 
Lift was obtained from the signal of electric strain gages mounted on 
the bearing supports. 

The power required to maintain a fixed amplitude of oscillation 
was measured visually with a wattmeter which Was incorporated in the 
electromagnetic— shaker circuit. At a corresponding time, signals from 
an ammeter and voltmeter in the same circuit were recorded on the 
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oscillogram. At the point of wing resonance in pitching, the product 
of the volts and the amperes was equal to the measured watts. 


METHOD 


The oscillating lifts were measured by strain— sensitive elements 
on the bearing supports. The incidence of the wing was indicated by 
the shearing strain in the torsion spring. The phase relation of the 
lift to the incidence of the wing was obtained from the simultaneous 
recording of the lift and angular position- on an oscillograph. The 
components of lift and the phase relation to wing incidence are 
illustrated in the following sketch. The frequency of oscillation is 
also obtained from the oscillogram. 



The aerodynamic damping moment, that is, the moment in phase with, 
the velocity, was obtained by the following procedures. Since a mome nt 
does work only with the component that is in phase with the angular 
velocity, the energy required to sustain an oscillation is a measure 
of the. damping moment and is given by rtQ 0 a 0 s ^- n 90* Therefore the 

input energy, obtained from the wattmeter and corrected for tare end 
electric circuit losses, yielded the damping moment. To confirm this 
value, the aerodynamic damping moment and tare moment were obtained at 
maximum response (sin cpg = 1) from the ammeter reading and the physical 

and electric constants of the machine. This aerodynamic moment could 
then be corrected for the tare damping moment. Two other methods of 
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measuring damping in common use were considered. If either the 
logarithmic decrement or the damping from the response curve is used, 
the aerodynamic damping moment can be obtained. The results of these 
methods were .in agreement. 

The aerodynamic moment in phase with the wing position may be 
obtained experimentally by observing the changes in the natural 
torsional frequency with Mach number. Such a method was used in 
reference 7 for 811 airfoil oscillating about its midchord. Since the 
model will oscillate in its '’damped" natural frequency when the forcing 
stimulus is removed, the inphase moment may be determined (appendix A). 
With this knowledge, the shift of the inphase center— of— pressure 
position from the quarter chord may be determined by dividing by the 
inphase lift . 


ACCURACY OF EXPERDffiNTAL RESULTS 


Before measuring the experimental oscillating lifts, the slope of 
the experimental static-lift curve of this wing was compared with the 
lift-curve slope for linear theory as demonstrated in figure 2. The 
theoretical dC^/da curve is the incompressible slope modified by the 


Glauert correction 



This corrected lift curve Is comparable 


to K' SI) for the case of = 0 as given by Dietze. Since there was 

some deviation of the experimental values from the theoretical values for 
the static case, the values of the experimental data for the oscillating 
case might also be expected to be less than the values indicated by the 
linear oscillatory theory. 


Quantitative measurements on high-frequency oscillations lead to 
certain problems not ordinarily encountered in static and low— frequency 
measurements. Deformation of the model under dynamic loading necessitates 
the correction of static calibrations. In the case of torsional 
oscillations the static calibration varied from the dynamic conditions 
as the result of the difference between the static and dynamic deflection 
curves. -This was corrected by applying a Holzer type of analysis to the 
dynamic configuration and relating the static strain relative to the 
midspan angle of incidence to the dynamic strain relative to the midspan 
angle of incidence. The deformation of the wing obtained from the Holzer 
analysis was approximately linear, the angle of attack of the side of the 
wing nearest the shakers being 9^.5 percent of the incidence at the 
opposite end for wing— rod 1. The angle of incidence used In the 
calculation of the aerodynamic coefficients is that at the midspan. 

In the measurement of the oscillatory lifts, the oscillograph 
records, particularly at the higher Mach numbers, were somewhat 
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distorted by wind— tunnel vibration and fluctuations in the flow of the 
air stream. To minimize the error due to this distortion, the oscillo- 
grams for both lift and angular position were read by taking an average 
over five consecutive cycles of an eight— point Fourier analysis. This 
analysis gave the amplitude of the lift vector and the magnitude of the 
angle of incidence as well as the phase between the lift and angular 
deflection. The precision of measurement as determined from repeated 
analysis from the same record gave lift values within 1 pound, the 
angle of incidence within 0.02°, and the phase angle within ±7°. Also, 
it was necessary to correct the phase relationship for any phase 
difference inherent in the galvanometer string. Unfortunately, the 
quadrature lift coefficient is very sensitive to variations of phase 
angle and therefore was subject to error. 

The correction due to dynamic magnification on wing— rod 1 was 
approximately 9 percent of the measured force (see appendix B). ho 
attempt was made to correct the experimental lifts and moments for the 
aerodynamic effect of this translation, since experiment and analysis 
showed the maximum displacement at midspan to be in the order of 
0.01 inch. This displacement for wing— rod 1 represents less than 
1 pound lift (approximately 1 percent total lift) under extreme tunnel 
conditions of high density and high Mach number. 

u 

There was a possibility of the damping moment being somewhat in 
error due to the variations of the tare damping. Since the magnitude 
of the tare was from 25 to 60 percent of the measured value, small 
variations in the tare led to large variations in the aerodynamic— 
moment coefficients. 

A small random error was introduced by the necessity of making 
interpolations and extrapolations of the theory in order to compare the 
theory with the experimental results. The interpolation was made by 
cross— plotting the coefficients against Mach number for a given value 
of reduced frequency. The extrapolations were made by extending the 
curves of coefficients against' reduced frequency given by Dietze and 
then cross— plotting the coefficients against Mach number at a given 
reduced frequency. This procedure tends to admit errors to the 
theoretical values. 


IMPROVEMENT OF TECHNIQUE 


For the future determination of oscillating air forces, there are 
various improvements to be effected. For more accurate results, 
dynamic calibrations for each specific test condition should be made. 
The difficulty in reading the oscillograms may be ameliorated by 
increasing the time scale or paper speed of the oscillograms. When 
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phase angles are being measured, it has been found that using the same 
type of signal in similarly damped galvanometer elements will reduce 
any phase difference due to the galvanometers. It is also important to 
check that the response curves of these elements are "flat" for the 
frequency range being investigated. 

The bending deflection of the wing— rod combination resulted in 
considerable dynamic magnification of the bearing reactions (see 
appendix B). The size of the correction is a function of the ratio of 
the forcing frequency to the natural bending frequency. To keep this 
magnification as low as possible, wings having very high bending 
frequencies compared to the forcing frequencies should be used. The 
tare damping in the system should be a minimum, the evaluation being 
made at all test amplitudes. 


DISCUSSION AND RESULTS 


The theoretical coefficients of lifts and moments were calculated 
by using the theory of Dietze (reference . l) . However, since the 
theoretical curves were determined to M = 0.7 only, all reference to 
theory above that Mach number refers to an extrapolation of these 
curves. For comparison of the lift coefficients it is well to note 
that the inphase and quadrature components are functions of both the 
magnitude and the phase of the lift vector. Under the test conditions 
it was impossible to treat these factors as independent variables. From 
figure 3 , it may be seen that a comparison can be made between the 
vector magnitudes of the experimental and theoretical lift coefficients. 
The difference between the experimental- and theoretical values is shown 
by the short vector marked "difference." It is shown that the 
experimental K'gp is insensitive to variations from the theoretical 

phase angle 9 -^. Conversely, the K"gp is highly sensitive to such 

variations. Since, the magnitude of the lift but not the phase angle 
could be determined accurately,- the experimental data which would be in 
closest agreement with theory would be magnitude of the total— lift 
vector. 

To show the range of values and agreement, the total theoretical 
lift has been plotted against experimental oscillating lift for wing- 
rod 1 (fig. 4). The ratio of experimental total lift to theoretical 
lift was plotted against Mach number (fig. 5)- 14 is clearly shown 

that the theoretical values are greater than the experimental values. 

It may be noted that the order of deviation from the theoretical values 
Is of the same magnitude as the deviation of the' static— lift data shown 
in figure 2. In the same graph is plotted 0 ^ showing the range of 

reduced frequencies in this set of tests. These and other data 
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pertaining to wing— rod 1, that is, Mach number, coefficients, and so 
forth, are also shown in table I. 

Since the theory assuming an infinite wing in an incompressible 
flow with some form of compressibility correction is in common usage, 
it is of interest to compare the experimental data with this theory. 

The ratio of experimental oscillating lift to this theoretical, lift 
plotted against Mach number is given in figure 6 . 

A comparison of the experimental phase angle with theory given in 
figure 7 shows that the trends indicated by theory are followed. When 
it is remembered that sin q>j_ is used in determining K" SI) (see fig. 3 ), 

it is only to be expected that the sign and order of magnitude of K" S p 
should deviate from the theory. Experimentally determined being 

insensitive to variations in the phase angle permits comparison with 
theory (fig. 8 ). It may be noted that the agreement is simiiar to that 
obtained for the magnitude of the lift vector. 

The opportunity was taken in this series of teats to examine the 
effects of small amplitudes of oscillation on the force coefficients. 

It has been suggested that for very small amplitudes some nonlinear 
effects associated with boundary layer may become significant. For the 
experiments herein, re ported no such effect was determined. The effect of 
various amplitudes of pitching oscillations is shown in figure 9 * It can 
be seen that the ratio of coefficients for the four Mach numbers shown is 
relatively constant and independent of amplitude, indicating that the 
amplitude of oscillation had little effect on the lift coefficient. 

In view of the large corrections in the lift, angle of incidence, 
and phase angle necessary in reudcing the data for wing— rod 3 , it was 
deemed advisable to omit the presentation of data and comparisons with 
theory. However, all trends were similar to those of wing— rod 1. 
Experiments were also conducted on wing— rod 2 which had a natural 
frequency of 4-5 cycles per second (midway between wing-rods 1 and 3 ). 

The data obtained for this combination (not presented) are subject to 
question. It was found that the measured value of damping in pitch 
decreased with Mach number, approaching' zero at M = 0 . 7 . The ratio 
of the experimental to the theoretical oscillation lift coefficients 
also approached a minimum at M = O. 7 . The data obtained for wing- 
rods 1 and 3 imply that this effect may be restricted to the particular 
test conditions for wing— rod 2. Unfortunately, time was not available 
to study this problem further,. It might be noted that the phenomenon 
exhibits certain characteristics of a self— excited vibration. The 
reason for the behavior is unexplained and it is hoped that future 
studies may answer the problem. 



10 


NACA EM L9K28a 


As previously mentioned, the experimental damping-moment coef- 
ficient in quadrature with the torsional displacement was determined. 

A comparison of these experimental results with theory plotted against 
Mach number (fig. 10) shows the experimental values to be consistently 
low. This result has also been noted by Bratt and Chinneck (refer- 
ence 7) for subsonic speeds. Some of the variation in the experimental 
coefficients may be attributed to fluctuations of the tare damping. 

The inphase moment coefficients K’pp were determined as shown in 

appendix A. Since, in the present case, the airfoil was oscillated 
about its quarter chord, the inphase moment was small. Thus, the 
change in natural torsional frequency may have been within the accuracy 
of the experiment. Therefore, a comparison for of experiment 

with theory shows little agreement (table I). However, for wing— rod 1 
which had the weakest spring and consequently the greatest sensitivity 
to aerodynamic moments, at the high Mach numbers (0.7 to 0.8) a 
negative oscillating moment occurred on the airfoil. This would indicate 
a rearward shift of the inphase center— of— pressure position of about 
3 percent chord. This center— of— pressure shift was predicted 
theoretically, but to a lesser degree (fig. 11). 


CONCLUDING REMARKS 


It is believed that the method of obtaining oscillating lifts and 
moments due to pitch as described in this paper can be refined to give 
sufficient accuracy for the experimental determination of those 
quantities in the transonic region. 

Good agreement of the experimental magnitude of the oscillating 
lift vector was obtained with the magnitude of the theoretical lift 
vectors of Dietze. Similar agreement was obtained between the 
experimental coefficients in phase with the wing displacement and the 
corresponding theoretical values. In general, both experimental 
coefficients were lower than the corresponding theoretical values. 

The trend with Mach number and reduced frequency indicated by 
theory for the phase between lift and incidence was obtained 
experimentally. The sensitivity of the quadrature components of lift 
to the variation in the experimental phase angle was such that the 
values obtained were questionable. 

The amplitude of torsional oscillation apparently has little effect 
on the lift coefficients in' the range of 0.004 to 0.04 radian. 
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A comparison of the experimental damping-moment coefficients with 
the theoretical values showed the experimental values to he consistently 
low. While quantitative agreement of the inphase moment coefficients 
with the theoretical values is not obtained, the movement of the inphase 
center of pressure tends to follow the behavior indicated by theory. 


Langley Aeronautical Laboratory 

Rational Advisory Committee for Aeronautics 
Langley Air Force Base, Va. 
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APPENDIX A 

DETERMINATION OF THE INPHASE MOMENT AND 
THE CENTER-OF-PRESSURE SHIFT 

The determination of the inphase moment, when the change of natural 
frequency of oscillation of a wing is known, is as follows: 

The definition of the moment is the air spring times the wing 
angular deflection, namely AKga 0 . The air spring is determined as 

follows: The natural frequency is 


% = 



the frequency with air spring force is 


“2 = 



+ ^9 
I 


thus 


then 



(Al) 


Now, since the air spring is known, the inphase moment is easily 
calculated. 

From inspection of equation (Al), when AKg ig positive, ojg > 
and there is a rearward shift of the center of pressure. When is 
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negative, > oig and there is a forward shift of the center of 

pressure. The position of the center of pressure referred to the 
quarter chord is determined by the ratio of the inphase moment to the 
inphase lift force, namely 


C.P. . .J*** 

ZLq COS <P;l 


This gives the center— of— pressure shift from quarter-chord position in 
percent chord. 

For comparison with theory, the center— of— pressure shift was 
calculated by the theory of Dietze (reference 3) as follows: 


C.P. 


Q 

Z'L 0 cos 

r2 


1 i s A - 

2 2 °\ 8 

pV 2 A*V 2 

2n l AZa 0 l — E 

2 °\ 2 


DD 


SD 


3 ,, 2 K* 

! 8^ “ K DD 


2 1 


2-r 


K’ 


SD 


O.375 
“’SD " ^ 
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DYNAMIC MAGNIFICATION OF LIFT FORCE 


When an elastic body is subjected to an oscillating force, it is 
possible to determine the amplitude response of the resulting vibration. 
In the case of a simply supported undam p ed beam carrying a uniformly 
distributed oscillating load (reference 8) it may be shown that 


where 



(Bl) 


and y s is approximately equal to the static midspan deflection under a 
uniform load W Q Z. In view of this relation, and recognizing that the 

deflection curve is a function of the loading, it follows that the 
reactions at the supports are also dependent upon the deflection curve. 
This fact may be shown as follows for a simply supported beam having an 
arbitrary loading: 



= = Ii +M -(eiJ 

^3 dx El dx dx 
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For constant El 


d 3 v = _1_ dM _ J_ 

3 El dx El 
dx 0 


where V, the vertical shear, is given hy E 
span loading. 



dx and \(x) is 


It may he shown that 


dim dV _ 

x-^0 ^3 ' El 


(B2) 


where E is the reaction at the support. Using the deflection curve 
from equation (Bl), the react ions, for a simply supported, undamped beam 
subjected to a uniformly distributed oscillating load are given by 


E 


= 


i= oo 

H 

i=0 


sin cot 


(.21 + l) 2 



(B3) 


Treatment of this nature roughly approximated the physical 
conditions of the experiment and gives an indication of the order of 
the correction. For the higher frequency combinations, a closer approach 
to actual conditions is obtained by considering the wing— rod as a two- 
span continuous beam built in at the wall. The reactions of this 
structure were determined by using the method of moment distribution. 

By assuming a given distributed load the deflection curve was obtained. 
The inertia effects w£y were then known for that loading. Using the 
combination of the distributed and inertia loads a new deflection curve- 
was obtained. By use of successive approximations the total loading was 
obtained. The ratio of the applied (aerodynamic) load to the measured 
load is then known. 

For the low— frequency wing— rod combination, this analysis, the 
results from equation '(B3) , and experimentally determined correction 
factors were in agreement.. The variation in the results of these 
methods for wing— rod 3j having a frequency ratio to/o^ of 0 . 65 , was 
such that only the order of magnitude of the correction factor could 
be determined. 
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Figure 1.— Diagrammatic view of wing in the 


Langley 4^ —foot flutter 
2 


research tunnel 







Figure 2.— Slope of the static— lift curve as a function of Mach number 

for wing— rod 1. 
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Figure 4.— Theoretical lift against experimental lift for wing— rod 1 
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Figure 5 .— Ratio of 
number 


experimental lift to theoretical lift against Mach 
and reduced frequency for wing— rod 1. 
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Figure 6.— Ratio of experimental lift to theoretical lift for 
incompressible flow at various Mach numbers and reduced 
frequencies. for wing— rod 1. 
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Figure 7»~ Theoretical and experimental phase-angle relationship to 
Mach number for wing— rod 1. (% = 1<?4 radians per second. 
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Figure 8.— Ratio of experimental K’g^ to theoretical K’gp against 
Mach number and reduced frequency for wing-rod 1. 
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Pitching displacement, radians 

Figure Ratio of experimental K'gp to theoretical K'g-p against 
pitching displacement for wing— rod 1. 
















